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Abstract 



' I We present a calculation of the 0{as) QCD corrections to the production of a Z boson in 

^^ 

Q ' association with two photons and to triple photon production at hadron colliders. All final-state 

f^, ' photons are taken as real. For the Z boson, we consider the decays both into charged leptons 

and into neutrinos including all off-shell effects. Numerical results are obtained via a Monte Carlo 

►^ , program based on the structure of the VBFNLO program package. This allows us to implement 

0> I general cuts and distributions of the final-state particles. We find that the NLO QCD corrections 

^J" ' are sizable and significantly exceed the expectations from a scale variation of the leading-order 

fT^ ■ result. In addition, differential distributions of important observables change considerably. The 

f~^ . prediction of two-photon-associated Z production with Z decays into neutrinos from the charged- 

^P , lepton rate works well, once we use an additional cut on the invariant mass of the charged-lepton 

pair. 

>'■ 
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I. INTRODUCTION 

The Large Hadron Collider (LHC) at CERN allows us to extend new-physics searches 
into energy regions never tested before. To claim a discovery of new phenomena, or provide 
improved bounds on them, it is important to know accurate predictions on their Standard 
Model (SM) background. This forces us to perform calculations beyond leading order, for 
integrated cross sections as well as for differential distributions. In this paper we present a 
calculation of the next-to-leading-order (NLO) QCD corrections to Z-boson production in 
association with two photons including off-shell effects for two decay modes 

.^ „ r "Z,77" pp,pp^Z77 + X->£+r77 + ^ (1-1) 

and to triple photon production 

"777" PP; PP -^ 111 + X . (1.3) 

Bothprocesses provide backgrounds for new-physics searches (for an overview, see e.g. 
Ref. pj). The process (II. 2p with its signature of two photons and missing transverse energy 
appears for example as background in models with gauge-mediated supersymmetry break- 
ing (GMSB) [2]]. In such models, the next-to-lightest supersymmetric particle is often a 
bino-like neutralino, which will decay into a gravitino and a photon. Since supersymmetric 
particles are pair-produced due to R-parity conservation, this leads exactly to the signature 
studied here. Triple-photonproduction provides a background to techni-pion production in 
association with a photon [3|, where the techni-pion decays into a photon pair. 

With this calculation the determination of triple vector-boson production cross sections 
at NLO QCD precision at hadron colliders is cornplete. The other processes belonging to 
this class have already been computed in Refs. |4J-|ll|. Multi- vector-boson production with 



one additional jet has been studied at NLO QCD so far for the diboson processes WWj, 
W'fj, WZj and ZZj jl2|-|lq and the triboson process W'yyj [ij 



The outline of the paper is as follows: in Section [IT] the Feynman diagrams which appear 
in our calculation are shown. We discuss the different techniques used to compute the virtual 
and real corrections and give an overview of the various checks which have been performed 
to ensure the correctness of the calculation. In section IIIII numerical results are presented. 
This includes scale variations of the leading-order (LO) and NLO integrated cross sections 
as well as selected differential distributions. We compare the decay of the Z boson into 
neutrinos, where the photons can only be radiated off the quark line, with the one into 
charged leptons, where also final-state photon radiation off the leptons is possible. A cut 
on the invariant mass of the lepton pair, restricting it to a small region around the Z mass, 
allows us to suppress the latter contributions. Section |IV] summarizes our work. 

II. CALCULATIONAL DETAILS 

We calculate all contributions up to order a^a^ to the processes (II. 11) and (II. 2p in the 
limit of massless fermions. In Fig. [H we show examples of Feynman diagrams appearing in 
our calculation. The different possibilities for the decay into charged leptons are depicted in 
the top row. Either all three vector bosons can be attached to the quark line (left), or one 
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FIG. 1: Examples of topologies of Feynman diagrams contributing to the different processes at 
tree level. Top row: pp — t- Z77 + X including Z decays into a pair of charged leptons. Bottom 
row: pp — )• Z77 + X including Z decays into a pair of neutrinos (left) and pp — )• 777 + X (right). 



(middle) or both (right) photons are radiated off the charged leptons. In total, 40 different 
Feynman diagrams appear at leading order. To improve th e sp eed of the calculation, we 
use the technique of "leptonic tensors" as described in Ref. [18|. Invariant subparts of the 
amplitude, namely the decays of the Z and virtual photons, are computed independently 
of the rest of the amplitude and only once per phase-space point. The result can then 
be easily reused when calculating the matrix elements, which greatly reduces the required 
computation time. Up to permutations of the vector bosons, only the diagram on the bottom 
left of Fig. [1] can appear at leading order for Z decays into a pair of neutrinos. Hence, the 
photons are always radiated off the quark line. 

The process ( II. 3p is correspondingly calculated up to order a^a^. Up to permutations of 
the photons, only the Feynman diagram depicted in the lower right of Fig. [1] contributes at 
leading order. 

To compute all these matrix elements, we use the helicity formalism of Ref. [13]. Contri- 
butions from virtual diagrams, which interfere with the tree-level diagrams, and from real 
emission occur at NLO QCD. Both are infrared divergent, but the sum of them does not con- 
tain any infrared divergences according to the Kinoshita-Lee-Nauenberg (KLN) theorem |20| . 
We achieve this cancellation in a numerically stable way by using the Catani-Seymour dipole 
subtraction method 2l| , which analytically extracts the divergences which are then canceled 
against the virtual ones. Part of the initial-state collinear singularities is factored into the 
parton-distribution functions, giving rise to additional so-called "finite collinear terms" . The 
NLO real corrections are obtained from the LO diagrams by attaching a gluon to the quark 
line in every possible way. This gluon can then be either a final-state particle, radiated off 
the quark line, or an initial-state one splitting into a quark-anti-quark pair. For the ^^77 
case there are 120 Feynman diagrams for each possibility. Hence, using leptonic tensors is 
particularly useful in this case. 

We obtain the virtual NLO QCD corrections by inserting gluon lines into the diagrams 
of Fig. [1] in every possible way. This yields loop diagrams with up to five external legs, 
i. e. pentagon diagrams. The occurring five-point integrals are solved using the prescription 



given in Ref. 



which does not need any inverse Gram determinants, with the setup of 



Ref. [23|. For loop integrals up to the box level, we use Passarino-Veltman reduction [24 
There we circumvent the explicit calculation of inverse Gram determinants numerically by 
solving a system of linear equations, which leads to a more stable behavior close to the 
critical regions. For all loop contributions, the full virtual corrections 
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factorize into a part proportional to the Born amplitude Mb and an infrared-finite remainder, 
My. Q denotes the partonic center-of-mass energy, i.e. the invariant mass of the photons 
and, for the Z'j'-f case, the leptons. We note that for diagrams where only a single vector 
boson is attached to the quark line (top right diagram of Fig. [1]), the virtual corrections 
completely factorize to their corresponding Born matrix element, i.e. they do not contribute 
to My. The finite remainders for all virtual corrections to two-boson amplitudes (top center 
diagram of Fig. [1]) are called "virtual-box" and to all three vector bosons attached to the 
quark line "virtual-pentagons" in the following. They were computed with the method 
described in Ref. |23|. To test the gauge invariance of the total Zj'y amplitude, we write 
the polarization vector of the Z as e^ = xzQz + ^z; where the last part simply denotes the 
remaining piece. Due to the structure of the amplitude, we cannot shift any contributions 
from pentagons to boxes, as has been done in other triboson processes |5|, la, Isl, |9[. Still, 
individual pentagon amplitudes change significantly depending on the value of xz- We have 
verified that their sum nevertheless stays unchanged. 

Isolated photons in the final state require particular care for the real emission part. 
Collinear emission of a photon from a massless quark leads to additional infrared diver- 
gences. Applying a naive separation cut between these two particles is not possible, how- 
ever, since this would restrict the phase-space of the final-state parton and therefore spoil 
the cancellation of divergences between virtual and real part. In this article, we solve the 
problem by imposing a particularly crafted cut 25|] between final-state parton (quark or 
gluon) i and photon •jj 



P% < PT-, 
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Rjii > ^0 ) 



(2.2) 



where j runs over all final-state photons and 6q is a fixed separation parameter. This cut 
allows the parton to be arbitrarily close to the photon, as long as its momentum vanishes 
simultaneously. This procedure avoids electroweak divergences, keeping at the same time 
the QCD pole of the real part intact. 

We have performed several cross checks to verify that our calculation is correct. First, 
all tree-level matrix elements have been compared against corresponding code generated 
by MadGraph |26|. We find an agreement at the level of the machine precision, i.e. at 
least 14 significant digits. Furthermore, we have compared the integrated cross sections of 
all leading-order processes as well as those with an additional jet against MadEvent and 
Sherpa [22] • For the triple photon process, internal technical cuts in these programs make 
the comparison difficult when we impose rather loose cuts on the photons. Therefore, we 
have performed an additional cross check with FeynArts 28|/FormCalc |29|/HadCalc 30 
for this instance. In all cases, we find an excellent agreement between the codes below the 
per mill level, compatible with the integration error. Also, the implementation of the Catani- 
Seymour subtraction scheme has been checked. We have tested that the large contributions 



from the true real-emission diagrams are canceled by the corresponding dipole terms once we 
approach the soft or coUinear regions. Additionally, we have verified that finite contributions 
appearing in the subtraction scheme can be shifted between the virtual and the real part 
without affecting the total result. Furthermore, for the virtual corrections we employ a 
gauge test for each phase-space point. 



III. NUMERICAL RESULTS 

We obtain numerical results from an NLO Monte Carlo program based on the structure 



of the VBFNLO program package [31] ■ In the electroweak sector, we take the masses of the 
W and Z boson as well as the Fermi constant as input. We then use tree-level relations to 
calculate the weak mixing angle and the electromagnetic coupling from these. As numerical 
values we hence use 



mw = 80.398 GeV 


sin^ {9w) = 0.22264 


niz = 91.1876 GeV 


a-^ = 132.3407 


6-^ = 1.16637 -10-^ GeV-2 . 





(3.1) 
We do not consider any effects arising from top quarks. All other quarks are considered to 



be massless. As parton distribution functions we use the CTEQ6L1 set at LO [32| and the 
CTIO set with as{mz) = 0.118 at NLO [33]. For both sets, we use five active flavors as 
initial-state quarks and in the renormalization group running of the strong coupling constant. 
As central value for the factorization and renormalization scale, we take the invariant mass 
of all uncolored particles, i.e. 

mz^^ = J{pe,/u^ +Ph/ur +P71 +^72)^ ^OT Z77 and 



m. 



777 



= V(P7i +P72 +^73)^ for 777- 



To take into account typical requirements of the experimental detectors, we impose the 
following set of minimal cuts on the transverse momentum and rapidity of the final-state 
charged leptons and photons 

Pt^ > 20 GeV pt, > 20 GeV \y^\ < 2.5 \ye\ < 2.5 . (3.3) 

Furthermore, photons, charged leptons and jets need to be well-separated in phase space, so 
they can be identified as separate objects in the detector. Also, divergences from coUinear 
photons need to be eliminated by cuts. Hence, we use additionally 

Ru > 0.3 R^^ > 0.4 Ri^ > 0.4 Rj^ > 0.4 Rj^ > 0.7 m^+^- > 15 GeV , (3.4) 

where the last requirement removes the singularity from a virtual photon splitting into a 
pair of oppositely charged leptons, when the invariant mass of the lepton pair becomes small. 
For the purpose of these cuts, a jet is defined as a final-state gluon or quark with px > 30 
GeV and \yj\ < 4.5. Furthermore, we impose the restriction defined in Eq. (12.21) with the 
parameter 60 = 0.7, which eliminates electroweak divergences between final-state quarks or 
gluons and photons coUinear to them. 

In Table [H we present results for the integrated cross section of Z77 and 777 production 
at the LHC with a center-of-mass energy of 14 TeV at the central scale hq. Besides our 



LHC 


LO [fb] 


NLO [fb] 


K-factor 


(t("Z77" -^ e+e~77) 








Pr^w > 20(20) GeV 


3.208 


4.877 


1.52 


PTj{e) > 30(20) GeV 


0.9987 


1.565 


1.57 


(t("Z77" -^ uv-yj) 








PT-y > 20 GeV 


2.905 


4.510 


1.55 


PTj > 30 GeV 


1.187 


1.906 


1.61 


cr("777") 








PTj > 20 GeV 


22.17 


58.57 


2.64 


PTj > 30 GeV 


6.907 


16.49 


2.39 



TABLE I: Total cross sections at the LHC for pp -^ Zjj + X with decays into charged leptons and 
neutrinos, and for pp — > 777. We show results for two sets of cuts at LO, NLO, and the associated 
K-factor. Relative statistical errors of the Monte Carlo are below 10"'^. 



Tevatron 


LO [fb] 




NLO [fb] 




K-factor 


(t("Z77" -^ 6^6^77) 








PT-y{e) > 10(10) GeV 


11.81 


16.37 


1.39 


PTj{e) > 20(10) GeV 


1.568 


2.226 


1.42 


(t("Z77" -^ uiy-yj) 








PT7 > 10 GeV 


2.421 


3.534 


1.46 


PT7 > 20 GeV 


0.5820 


0.8492 


1.46 


a("777") 








PT-y > 10 GeV 


40.67 


84.20 


2.07 


PT7 > 20 GeV 


6.008 


10.28 


1.71 



TABLE II: Total cross sections at the Tevatron for pp —?■ Z77 + X with decays into charged 
leptons and neutrinos, and for pp — )• 777. We show results for two sets of cuts at LO, NLO, and 
the associated K-factor. Relative statistical errors of the Monte Carlo are below 10~^. 



standard set of cuts in Eqs. 03. 3p and (13. 4p . we also show results with a modified cut on 
the transverse momentum of the photon of 30 GeV. The results for the Tevatron with its 
center-of-mass energy of 1.96 TeV are denoted in Table HTl Here, we have reduced the cut on 
Pti to 10 GeV, and besides the standard value of 20 GeV also give results for px^ larger than 
10 GeV. In all cases, we consider only a single generation of charged leptons or neutrinos in 
the final state. 

From hereon, we will restrict ourselves to results for the LHC. Additionally, we now sum 
over both electrons and muons in ^^77 and over all three neutrino generations in .^1,77. 
In Figs. |5] and [31 we show the dependence of the cross section on the factorization and 
renormalization scale for the ^'77 process with decays into charged leptons and neutrinos, 
respectively. We vary both scales together, independently and inversely, such that when one 
scale is increased, the other one decreases by the same factor, in the range 
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FIG. 2: Left: Scale dependence of the total LHC cross section for pp — > ^77 + X ^ £"'"^"77 + X 
at LO and NLO within the cuts of Eqs. (|3.3l 13. 4p . The factorization and renormalization scales 
are together, inversely or independently varied in the range from 0.1 ■ //q to 10 • /Uq- We show NLO 
curves without and including an additional veto on -pTj of 50 GeV. Right: Same as in the left panel 
without jet veto, but for the different NLO contributions at fip = fJ-R = Cmo with /^q = mz-yy- 
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FIG. 3: Left: Scale dependence of the total LHC cross section for pp — ?> .Z'77 + X ^- 2/^77 + X 
at LO and NLO within the cuts of Eqs. ()3.3l 13. 4p . The factorization and renormalization scales 
are together, inversely or independently varied in the range from 0.1 ■ /io to 10 • //q- We show NLO 
curves without and including an additional veto on pTj of 50 GeV. Right: Same as in the left panel 
without jet veto, but for the different NLO contributions at fip = /i/j = ^/Uq with /xq = rnz-y-y- 



where the central scale /xq has been given in Eq. (13.21) . The corresponding plot for the 777 
process is displayed in Fig. |H In all cases, we also show results where we have applied an 
additional veto on jets with pxj > 50 GeV. This results in smaller renormalization scale 
variations, but the principal shapes of the different curves are not affected. We observe that 



80 
70 
60 
50 



=- 40 



30 
20 
10 



..^ rn 


y' 

,y' 
■■■■.., NLO ^^^ 






y^' , ' ' 


NLOp\,„,= 50GeV - 




^ 


^"'^ 


solid: Uf = I^R = 5 M^ 
dashed:np = 5M^,Hpi=M^ 
dotted: Up = M^, jif,= ^M^ 
dotted-dashed: Hp = ^ M^, np= 1/^M^ 



0.1 



10 



5 





yyy 




^P = Hr = 5 M^ 


70 


" 




Total NLO 


60 






50 


- 




- 


£40 


^.,^ 




_ 








Virtual-born 


30 




Real 


...;:..:^--:..-;.;;"' 


20 





'"- 


10 


- 




- 








Virtual-Pentagons 



0.1 



10 



5 



FIG. 4: Le/it; Scale dependence of the total LHC cross section for pp -^ 777 + X at LO and 
NLO within the cuts of Eqs. (J3.31 13. 4p . The factorization and renornialization scales are together, 
inversely or independently varied in the range from 0.1 -/xq to 10 -/io- We show NLO curves without 
and including an additional veto on pj-,- of 50 GeV. Right: Same as in the left panel without jet 
veto, but for the different NLO contributions at fip = /^K = S.fJ-o with /xq = m^^^. 



for all processes, the NLO corrections are much larger than estimated by a scale variation 
of the LO result. The respective K-factors at the central scale are 1.52 and 1.55 for ^'77 
decaying into charged leptons and into neutrinos, and 2.64 for 777. We observe that the 
NLO cross section typically rises with the factorization scale and decreases with the renor- 
malization scale. This leads to cancellations for the combined scale dependence which can 
yield a very fiat behavior, as can be seen for example in Fig. [2l In contrast, when we vary 
the two scales inversely, the variations add up and lead to a rather significant dependence. 
This almost fiat behavior is therefore accidental, and should not be considered as a very 
small remaining scale uncertainty. 

On the right-hand side of Figs. [2], [3] and HI we show the individual contributions to 
the unvetoed NLO cross section for the different processes as function of the combined 
factorization and renormalization scale. Both the real-emission part, which contains the true 
real-emission contribution, the dipole terms from the Catani-Seymour subtraction scheme 
and the finite coUinear terms, as well as the virtual part proportional to the Born matrix 
element exhibit a strong scale dependence, which partly cancels in the sum. For the 2^77 
process, the bulk of the contribution is given by the virtual part, while for 777 it is evenly 
distributed between the two. This is due to the higher average partonic center-of-mass energy 
for ^77, which reduces the initial-state gluon parts due to the steep fall-off of the gluon pdfs 
with growing energy. The pentagon finite virtual remainders, "virtual-pentagons" , defined 
in Eq. (12. ip . are small in all cases. Finite box virtual remainders, "virtual-box", exist only 
for the Zi'-ff process and they are small as well. 

Finally, we look at several differential distributions for all considered processes. In each 
figure, we show at the left-hand side the differential cross section both at LO and NLO. The 
latter we plot without and including an additional veto on jets with ptj > 50 GeV. On 
the right-hand side, we depict the differential K-factor for both the unvetoed results and 
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FIG. 5: Left: Transverse- momentum distribution of the photon with largest transverse momentum 
in Zi'j^ production with Z decaying into charged leptons for the LHC. We show LO (dotted blue 
line) and NLO cross sections without (solid red) and including (dashed black) a jet veto of 50 GeV, 
using the cuts of Eqs. (13.31 l3-4p . Right: Associated K-factor as defined in Eq. ()3.6p without (solid 
red) and including (dashed black) the jet veto. 
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FIG. 6: Left: Transverse- momentum distribution of the photon with largest transverse momentum 
in Ziyjj production with Z decaying into neutrinos for the LHC. We show LO (dotted blue line) 
and NLO cross sections without (solid red) and including (dashed black) a jet veto of 50 GeV, 
using the cuts of Eqs. (|3.3l 13. 4p . Right: Associated K-factor as defined in Eq. (|3.6p without (solid 
red) and including (dashed black) the jet veto. 



including the jet veto. It is defined as 
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FIG. 7: Left: Transverse- momentum distribution of the photon with largest transverse momentum 
in 777 production for the LHC. We show LO (dotted blue line) and NLO cross sections without 
(solid red) and including (dashed black) a jet veto of 50 GeV, using the cuts of Eqs. (J3.31 13. 4p . 
Right: Associated K-factor as defined in Eq. (l3.6p without (solid red) and including (dashed black) 
the jet veto. 
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FIG. 8: Left: Transverse- momentum distribution of the photon with largest transverse momentum 
in .^^77 production with Z decaying into charged leptons for the LHC. We show LO (dotted blue 
line) and NLO cross sections without (solid red) and including (dashed black) a jet veto of 50 GeV, 
using the cuts of Eqs. ()3.3l 13. 4p and the additional cut on mu given in Eq. (j3.7p . Right: Associated 
K-factor as defined in Eq. (|3.6p without (solid red) and including (dashed black) the jet veto. 



where x denotes the considered observable. 

In Figs. [5], [6] and [TJ we present the transverse-momentum distributions of the photon 
with the largest transverse momentum for Zi'-ff, Z^^^ and 777, respectively. For each of 
the different processes, the K-factor for the differential distribution is not constant. Rather, 
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FIG. 9: Left: Distribution of missing transverse momentum in Z^^'j production with Z decaying 
into neutrinos for the LHC. We show LO (dotted blue hne) and NLO cross sections without (soHd 
red) and including (dashed black) a jet veto of 50 GeV, using the cuts of Eqs. (J3.31 I3.4p . Right: 
Associated K-factor as defined in Eq. (|3.6p without (solid red) and including (dashed black) the jet 
veto. 
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FIG. 10: Left: Separation of the two photons in Zy^^ production with Z decaying into neutrinos 
for the LHC. We show LO (dotted blue line) and NLO cross sections without (solid red) and 
including (dashed black) a jet veto of 50 GeV, using the cuts of Eqs. (|3.3ll3.4p . Right: Associated 
K-factor as defined in Eq. (|3.6p without (solid red) and including (dashed black) the jet veto. 

we observe for the Z'j'j processes a significant rise with larger transverse momenta. At small 
transverse momenta, the K-factor is close to the integrated results, as this is where the bulk 
of the cross section lies. Once we apply the additional jet veto, the K-factor instead drops 
to unity as the transverse momentum of the photon increases. Hence, the large differential 
K-factors are caused by events where the leptonic system recoils against the additional jet. 
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For tri-photon production (Fig. [7]), we observe a substantially larger K-factor for very small 
transverse photon momenta, where the recoil against the jet helps to fulfill the transverse- 
momentum cut, but is not yet large enough to trigger the possible jet veto. Once we go to 
large values, the K-factor is almost constant without jet veto, and gradually drops to one 
including it. 

Furthermore, we show this distribution in Fig. [H] for ^£77, where we have applied the 
additional cut 

mz - ITz < mu < mz + 2Tz (3.7) 

on the invariant mass mu of the lepton pair. This cut allows us to restrict the process in 
such a way, that the photons are predominantly radiated off the quark lines (top left diagram 
of Fig. [1]), which is the only possibility for Z decays into neutrinos. Comparing this figure 
with Fig. ini we see that the shape of the plots both for the cross sections as well as for 
the K-factors is almost identical. We have also checked several other distributions. Up to 
a global normalization constant, there are only small and negligible differences between the 
two. Therefore, by using this cut it is possible to extrapolate the differential rate with Z 
decays into neutrinos from the charged-lepton rate. 

In Figs. [9] and [TOl we present two more distributions for the ^'77 process with decay into 
neutrinos. The distribution of missing transverse momentum originating from the neutrino 
pair is shown in Fig. |9l At LO, we see a broad distribution superimposed with a peak above 
2pT^,min = 40 GeV, which is smeared out at NLO. This gives rise to the structure at low 
momenta in the K-factor, which we observe in the figure. For large missing momentum, we 
see again a significantly larger K-factor without any veto. With an additional jet veto of 50 
GeV the K-factor drops to unity in this region. 

Finally, in Fig. [TU]we show the separation in the rapidity-azimuthal-angle plane between 
the two photons, i?^-y, for ^1,77 . Also here, we observe a large variation of the differential 
K-factor with the value of the i?-y-y separation. Both curves without and including the 
additional jet veto exhibit a similar shape. 

Therefore, it is necessary to calculate fully differential NLO cross sections to obtain the 
correct shape of the distributions. A simple rescaling of the LO result with the integrated 
K-factor would give large deviations. 



IV. CONCLUSIONS 

We have calculated the NLO QCD corrections to the processes pp, pp — )■ Z77 + X with 
decays of the Z boson into both charged leptons and neutrinos and including all off-shell 
effects, as well as to pp,pp — )■ 777 + X. These processes appear as important background 
in searches for new physics, for example in supersymmetric GMSB models or the latter in 
technicolor models. 

In all cases, we obtain large corrections to the leading-order cross section. The typical 
size of the integrated K-factors is around 1.5 for Z'-f-f and 2.6 for triple-photon production. 
This strongly exceeds the expectations from LO scale variation, which is at the level of 10% 
to 16%. For an estimate of the remaining scale dependence, we vary the factorization and 
renormalization scale by a factor two around the central scale /xq = ''^277/777 such that the 
product of the two scales stays constant, i.e. when one scale takes the factor two, the other 
has factor | and vice versa. This gives the largest differences between ^ = 2 and ^ = | with 
changes of 9.0%, 7.5% and 18.6% for ^^77, ^',^77 and 777 at NLO, respectively. For a joint 
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variation, the difference is below 4% for all processes due to accidental cancellations. 

The size of the NLO corrections exhibits a strong dependence on the observable and on 
different regions of phase-space. We also observe that for the processes considered in this 
article an additional jet veto does not flatten the differential K-factors. Hence, it is necessary 
to have a dedicated fully-exclusive NLO parton Monte Carlo program available. We have 
also studied, whether a measurement of the charged-lepton decay mode in ^'77 will allow 
us to predict the neutrino case. Here, we find that this works well with an additional cut, 
which restricts the invariant mass of the lepton pair to the region close to the Z pole. 

With these processes, we have completed the NLO QCD calculations of triple vector- 
boson production at hadron colliders including off-shell effects, Higgs resonances and leptonic 
decays. All these processes will be available in the next release of VBFNLO. 
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